Controllable soliton emission from a Bose-Einstein condensate. 
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We demonstrate, through numerical simulations, the controllable emission of matter-wave bursts 
from a Bose-Einstein Condensate in a shallow optical dipole trap. The process is triggered by 
spatial variations of the scattering length along the trapping axis. In our approach, the outcoupling 
mechanism are atom-atom interactions and thus, the trap remains unaltered. Once emitted, the 
matter wave forms a robust soliton. We calculate analytically the parameters for the experimental 
implementation of this atomic soliton machine gun. 
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PACS numbers: 42.65.Jx, 42.65.Tg 

Introduction.- After the remarkable experimental real- 
ization of Bose-Einstein condensates (BEC) with alkali 
atoms [l J it was soon realized that the coherent behavior 
of the atom cloud could be used as the basis for an "atom 
laser". The first realization of such device used short 
radio-frequency pulses as an outcoupling mechanism, 
flipping the spins of some of the atoms and releasing them 
from the trap |2|. Later, other atom lasers were built 
with different configurations leading to pulsed, semi- 
continuous or single-atom coherent sources^ [3. 

There are important differences between these devices 
and photon pulsed lasers. The nature of the waves and 
the lack of a population inversion mechanism are the 
most evident. However, from the practical point of view, 
one of the most remarkable is that a pulsed BEC laser 
will show a significant spreading of the atom cloud if the 
number of atoms per pulse is large. 

A way to overcome this difficulty would be to use 
atomic bright solitons as matter wave pulses. This 

idea was explored in Ref. [10( to generate a train of soli- 
tons by the mechanism of modulational instability (Ml) 
and was referred to as an atomic soliton laser. However, 
although one could extract a few coherent solitons from 
such a system this atomic soliton laser would be very 
limited since: (i) the final number of atoms per pulse af- 
ter MI is only a small fraction of the initial number of 
atoms in the condensate due to collapse processes, (ii) 
the number of solitons generated is not large and half of 
them would be directed backward (iii) the trap should be 
destroyed for outcoupling and (iv) the pulses will travel 
at different speeds once the trap is removed. 

Thus, it is important to discuss new outcoupling mech- 
anisms for atom lasers. This is specially interesting since 
the techniques for generating BEC with growing number 
of particles and their physical properties are nowadays 
well established and the current challenges in the field 
face the applications of coherent matter waves to the de- 
sign of practical devices 



In this letter, we show how a highly controllable train 
of up to several hundreds of matter-wave solitons can be 
extracted from a BEC without altering the trap properties 
but instead acting on the scattering length, by changing 
it along the atom cloud. 

System configuration and theoretical model.- Let us as- 
sume that a large BEC is strongly trapped in the trans- 
verse directions (x, y) and weakly confined in the longitu- 
dinal one (z). We now consider the effect of a sharp vari- 
ation along z of the scattering length, which is changed 
from positive (or zero) to negative, thus making it in- 
homogeneous. This can be done with magnetic |l2| or 
optical^ techniques and therefore this region of nega- 
tive scattering length can be displaced along the conden- 
sate. Thus, when it is located close to one edge of the 
trap, overlapping the wing of the cloud, the tail of the 
condensate may be able to form a single soliton which, 
because of its higher internal energy will be outcoupled 
from the cloud and thus emitted outward. When the 
condensate refills the gap left out by the outgoing pulse 
a new soliton would be emitted. This process would con- 
tinue while there is a large enough remnant of atoms in 
the trap and would lead to a soliton burst escaping from 
the BEC. 

Thus, we will describe our system of N weakly inter- 
acting bosons of mass m, trapped in a potential V(r) is 
the mean field limit by a Gross-Pitaevskii equation of the 
form 



lh ~dt 



h 2 

-— V 2 * + y(f)* + C/(z)|*| 2 *, (1) 
2m 



where ^ is the condensate wavefunction, N = J |\]/| 2 d 3 r 
The coefficient U(z) — Anti 2 a(z) / 'm characterizes the 2- 
body interaction and since we will consider spatially in- 
homogeneous systems it will be a function of z. 

In this paper, we consider a BEC tightly confined in 
(x, y) by a harmonic potential V± and more relaxed along 
z where there is and optical dipole trap V z [3E3. Thus, 



2 



we have 



V(r) = V x +V z 



(x 2 + z/ 2 )+y 



1 — exp 



z 

L2 



(2) 

where Vn is the depth of the shallow optical dipole po- 
tential and L is the characteristic width of the trap. We 
will consider situations in which the ground state of the 
optical dipole trap is much larger than the ground state 
of the transverse harmonic potential. 

In FigQ] we have plotted the geometry of the system 
The choice of a shallow Gaussian trap is very impor- 
tant for our model, since we are interested on studying 
the outcoupling of solitons along the z axis. Thus, we 
need a potential barrier that can be overcome by the 
self-interaction effects. In this situation we can describe 
the dynamics of the condensate in the quasi-one dimen- 
sional limit as given by a factorized wavefunction of the 
form 0| ^(r,t) = 3>a(x,y) ■ tp(z,t), satisfying 



«^ = -^0+/W + ^lV, (3) 



where rj_ = y/ K/rrwx is the transverse size of the cloud, 
f{rj) — V z /(hiy±/2), t = v±t is the time measured in 
units of the inverse of the radial trapping frequency and 
g{z) = ^/8Trr^_a(z) is the effective interaction coefficient. 




FIG. 1: Sketch of the system showing the axial trap (dashed 
line) and the shaded region z > D (=2L in this particular 
case) which has negative scattering length a. For the rest of 
the cloud (z < D) a is set to zero. The black zone displays 
the overlapping of the cloud and the region with negative a. 
This causes the outcoupling for a critical value of the scat- 
tering length a cr . The continuous lines are the numerically 
calculated profiles of eigenstates for a) a» = , b) 
and c) a» = 1.2a cr . 

In Fig. ^we show a sketch of the setup to be considered 
in this paper. We show the axial optical dipole potential 
of depth Vq and half width L. The shaded region (z > 



D) has negative scattering length. A spatial variation 
along z of the scattering length a, can be achieved by 
magnetically tunning the Feschbach resonances or 
by their optical manipulation by means of an additional 
laser beam . In this paper we choose a step model for 
a(z) in order to present the basic mechanism. We must 
stress that our ideas apply to smoother variations of a. 
Thus, we will use a dependence of the form 



a(z) 



0, z < D 

a* < 0, z > D. 



(4) 



The continuous lines in Fig. Q]are the numerically calcu- 
lated profiles of eigenstates of Eq. © for three different 
values of a*. The lower curve corresponds to he linear 
case (a* = 0), whereas lines b) and c) display the shape 
of the cloud for and a* = 1.2a cr , respectively. 

The value CL(2.f IS the (negative) scattering length needed 
for emitting one soliton from the cloud. We will calculate 
it by means of an approximate analytical technique. 

Single soliton emission.- In Fig. 1 (black region), it is 
shown the overlapping of the tail of the trapped cloud 
and the zone where a = a* . Due to the atom-atom inter- 
action, the trapping changes in this region. Thus, if a* is 
negative enough, the minimum of the effective potential 
will be displaced toward this region and the cloud will 
move to the shaded region of Fig. ^ 
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FIG. 2: [Color Online] Effective potentials corresponding to 
different values of the scattering length: (a) a, — 0, (b) a* = 
a cr , (c) a„ = 1.2a cr . In (d) we show a pseudocolor plot of 
|i/j(z,t)| 2 showing the soliton emission for a» = 1.2a cr . The 
parameter values are Vb = hv±/2, v± = IK Hz, L — 4r±, 
D — 2L, w = 5Ar± and N ■ a cr — — 100/im. The spatial 
window size is 100r^, and time goes from to 4000I/J 1 . 

To study this possibility we will use an averaged La- 
grangian formalism |l7| for studying localized solutions 
of Eq. © . We will use a Gaussian ansatz 



il){z,T)=Ae-( z - z ^l 2w \ i ^ z \ 



(5) 
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where z q (t) accounts for the motion of the center the 
cloud with velocity proportional to v. The standard cal- 
culations lead to the equations 



n(*o) 




(6a) 



(6b) 



where erfc(u) = f u exp(— v 2 )dv is the complemen- 
tary error function. Thus, zq evolves like a classical par- 
ticle under the effect of a potential II(zo)-This provides 
a qualitative understanding of the soliton emission: for 
the linear case (a = 0) the center of the cloud is located 
at the bottom of the Gaussian trap, which corresponds 
to its fundamental eigenstate [Fig. ^ (a)]> as a* takes 
more negative values the effective trapping of the cloud 
is deformed and the minimum of the equivalent potential 
moves to the region with z > 0. Fig. [5] shows the equiv- 
alent potentials II given by Eq. I|6bjl for different values 
of a* < 0. Fig. Efa) shows the effective potential for 
a* = 0. As a* is decreased there is a limiting value a cr 
for which the potential 11(0) = IT(oo) [Fig. Ejfb)], thus if 
the atom cloud is initially placed at zq — it will oscillate 
around the minimum and escape zq(t) — > oo for r — > co, 
a phenomenon which is called soliton emission^}. The 
critical value of a* that corresponds to the threshold for 
soliton emission can be obtained within our formalism 
from the condition, 11(0) = II(oo) which yields to: 



Na r _ 



8nVo wL 



erfc 



\/2D 



wL 



hv± \/L 2 + w 2 



L 

5. 0A-, 

w 



(7a) 



(7b) 



the first approximation is valid for D > w and the second 
for Vo = Ahv±, with A < 1, which provides a shallow 
trap with w L. We have found that the ratio between 
the exact numerical value and the prediction from Eq. 
(|7a|l is a factor 1.15 which is essentially due to the specific 
choice of the ansatz. In fact, once the soliton is emitted 
the wavefunction takes a hyperbolic secant shape and 
its width differs from the Gaussian by a similar factor. 
However, taking a sech profile as variational ansatz does 
not yield to analytical results with the potentials we have 
considered. 

For more negative values of a < a cr not only the soliton 
is outcoupled from the system but also it propagates with 
a finite asymptotic speed in the nonlinear medium. An 
example of the effective potential for a — 1.2a cr is shown 
in Fig. |2fc) and the corresponding soliton emission in 
Fig. [2d) , for a 7 Li condensate such as that of Ref. Q 



Partial outcoupling and multisoliton emission.- A 
deeper numerical exploration based on Eq. reveals 
more interesting effects which are beyond the averaged 
lagrangian description. An example is shown in Fig. 3 
for parameter values Vo = Hu±/2, vy_ = IK Hz, L = Ar±, 
D = 2.5L, N = 3 • 10 5 , w = 5.4r x . The vertical axis in 
each figure is time from t = to t = Is and the hori- 
zontal width of each window is 100r^ 300/xto. When 
a* = 0.9a cr the atom cloud is only slightly deviated to 
the region with a < [Fig. El (a) ] . Decreasing a* down 
to a* = 1.95a cr leads to the phenomenon of soliton emis- 
sion with some remnant staying on the optical dipole 
trap [Fig. Efb)]. However, as a* is made more and more 
negative we obtain the emission of an integer number of 
solitons. 
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FIG. 3: [Color Online] Multiple emission of atomic solitons 
from a BEC reservoir for different values of a,: (a) a, = 
0.9a cr , (b) a„ = 2.0a cr , (c) a„ = 3.3a cr , (d) a* = 4.8a cr , (e) 
a* = 6.5a cr , (f) a* = 8.7a cr .The vertical axis in each figure 
is time from t — to t — Is and the horizontal width of each 
window is 100rx ~ 300/im. The rest of the parameters are 
the same as for the inset of Fig. |21 except D = 2.5L. 
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FIG. 4: Soliton emission for D = 2.5L and five different times 
from to Is. Left: Emission of a single soliton for a = l.la cr - 
Right (b) Emission of five solitons for a — 4Aa cr - 



The phenomenon of partial emission can be qualita- 
tively understood with the variational method. As it can 
be seen in Fig. Elthe effective potentials n have only one 
minimum for D = 2.0L. However if the value of D is in- 
creased to say D = 2.5L, the potential n for values close 



4 



to a cr presents one maximum between two adjacent min- 
ima. This configuration of the effective potential causes 
the split of the matter wave in two parts: one remains in 
the trap and the other is emitted as a soliton. 
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FIG. 5: Average time interval r so ; between the emission of 
two solitons as a function of a,/a cr . The curve follows a law 
of the type r oc a~ 2 (solid curve). Insert: Number of solitons 
emitted as a function of a»/a cr - 

From our numerical exploration we have found some 
regularities that are noticeable. In first place, as it can 
be appreciated from Fig. the velocities of the outgo- 
ing pulses are almost equal. Secondly, the intervals of 
emission are very regular and decrease with the effec- 
tive nonlinearity with a law of the form r ~ (a*/a cr )~ 2 . 
Finally, the values of the scattering length required for 
outcoupling several solitons are almost exact multiples of 
a cr N as it is clear by the approximate linear dependence 
of the number of solitons emitted in the inset of Fig. 03 
This can be interpreted as the requirement of a critical 
effective nonlinearity remaining in the trapped cloud for 
the formation of every soliton. Finally, we must comment 
from Fig. 0] that the amplitude of the emitted solitons 
decays from the first one of the burst to the last approx- 
imately in a 30%, depending on the number of pulses. 

Discussion and conclusions.- These regular and con- 
trollable soliton trains could be realized with any atomic 
specie for which the scattering length could be made neg- 
ative on a region of the space. This has been done with 
several atomic species such as 7 Li, 85 Rb and 133 Cs. Tak- 
ing for instance the experimental parameters of Rcf. 8] 
for Lithium and assuming 3 x 10 5 atoms in the optical 
dipole trap we obtain that a few (2-3) solitons could be 
emitted for a depth of Vq = Hv±/2. However, increasing 
the initial particle number in the zero-scattering length 
region could allow the observation of soliton trains with 
a higher number of solitons, e.g. for 3 x 10 6 we estimate 
up to N so i ~ 20. For 85 Rb and taking parameter data 



from Ref. we estimate about 25. Even more inter- 
esting effects could be obtained with Cesium, for which 
experimental data of Ref. [13] lead to an estimate of 
N so i ~ 200. In all cases these numbers can be raised or 
lowered by acting on the parameters: a*, N, D and Vq. 

In summary, we have proposed a novel mechanism for 
outcoupling coherent matter wave pulses from a Bosc- 
Einstein Condensate. Our system is able to perform a 
regular and controllable emission of atomic soliton bursts 
that are easily extracted by an adequate choice of the 
control parameters. Using this mechanism a train of even 
several hundred of solitons could be coherently outcou- 
pled from a condensate. As the techniques for coher- 
ently feeding the remaining condensate progress our idea 
could provide an outcoupling mechanism for a continuous 
atomic soliton laser. 
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